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Abstract. For any t > 0, we present an algorithm which takes as input 
a semi-algebraic set, S, defined by Pi < 0,...,P S < 0, where each Pi € 
R[Xi, . . . has degree < 2, and computes the top I Betti numbers of S, 

fefc_l(S), . . . , bk—e(S), in polynomial time. The complexity of the algorithm, 
stated more precisely, is 5Zi=o (i)^ 2 °* . For fixed £, the complexity 

of the algorithm can be expressed as s^+ 2 fc 2 °' ' , which is polynomial in the 
input parameters s and k. To our knowledge this is the first polynomial time 
algorithm for computing non-trivial topological invariants of semi-algebraic 
sets in R fe defined by polynomial inequalities, where the number of inequalities 
is not fixed and the polynomials are allowed to have degree greater than one. 
For fixed s, we obtain by letting £ = k, an algorithm for computing all the 
Betti numbers of S whose complexity is k 2 



1. Introduction 

Let R be a real closed field and ScR'a semi-algebraic set defined by a Boolean 
formula with atoms of the form P > 0, P < 0, P = for P e V C R[X U ...,X k ]. 
It is known H2 HH HE El [E] that the topological complexity of S (measured 
by the various Betti numbers of S) is bounded by 0(s 2 d) k , where s = #('P) and 
d = maxpgp dcg(P). Note that these bounds arc singly exponential in k. More pre- 
cise bounds on the individual Betti numbers of S appear in [5] . Designing efficient 
algorithms for computing the homology groups, and in particular the Betti num- 
bers, of semi-algebraic sets are considered amongst the most important problems 
in algorithmic semi-algebraic geometry. 

Even though the Betti numbers of S are bounded singly exponentially in fc, 
there is no known algorithm for producing a singly exponential sized triangulation 
of S (which would immediately imply a singly exponential algorithm for comput- 
ing the Betti numbers of S). In fact, the existence of a singly exponential sized 
triangulation, is considered to be a major open question in real algebraic geometry. 
Doubly exponential algorithms (with complexity (sd) 2 ' k ) for computing all the 
Betti numbers are known, since it is possible to obtain a triangulation of S in dou- 
bly exponential time using cylindrical algebraic decomposition. In the absence of 
singly exponential algorithms for computing triangulations of semi-algebraic sets, 
singly exponential algorithms are known only for the problems of testing emptiness 



Key words and phrases. Betti numbers, Quadratic inequalities, Semi-algebraic sets. 

The author was supported in part by an NSF Career Award 0133597 and a Sloan Foundation 
Fellowship. A preliminary version of this paper appears in the Proceedings of the ACM Symposium 
on Theory of Computing, 2005. 



2 



SAUGATA BASU 



[27(17]. computing the zero-th Betti number (i.e. the number of semi-algebraically 
connected components of S) [HI [2Ql El HH [8] , as well as the Euler-Poincare char- 
acteristic of 5* [1] . Very recently a singly exponential time algorithm has been given 
for computing the first few Bctti numbers of semi-algebraic sets [3] (see also [6]). 

In this paper, we consider a restricted class of semi-algebraic sets - namely, 
semi-algebraic sets defined by a conjunction of quadratic inequalities. Since sets 
defined by linear inequalities have no interesting topology, sets defined by quadratic 
inequalities can be considered to be the simplest class of semi-algebraic sets which 
can have non-trivial topology. Such sets are in fact quite general, since every 
semi-algebraic set can be defined by (quantified) formulas involving only quadratic 
polynomials (at the cost of increasing the number of variables and the size of the 
formula). Moreover, as in the case of general semi-algebraic sets, the Betti numbers 
of such sets can be exponentially large. For example, the set ScR defined by 

X 1 (l-X 1 ) <0,...,X k (l-X k ) <0, 

has b (S) = 2 k . 

Hence, it is somewhat surprising that for any fixed constant £, the Betti num- 
bers bk-i(S), . . . ,bk-e{S), of a basic closed semi-algebraic set S c R h defined by 
quadratic inequalities, arc polynomially bounded. The following theorem appears 
in0. 

Theorem 1.1. Let R a real closed field and S C R fe be defined by 
Pi < 0, . . . , P 3 < 0, dcg(P) < 2, 1 < i < s. 

Then, for £>Q, 

6*^(5) <(J)*° W - 

Notice that for fixed £ this gives a polynomial bound on the highest I Betti 
numbers of S (which could possibly be non-zero). Observe also that similar bounds 
do not hold for sets defined by polynomials of degree greater than two. For instance, 
the set VcR* defined by the single quartic inequality, 

k 

J^XfiX,- l) 2 -e>0, 

i=l 

will have b k -i(V) = 2 fc , for all small enough e > 0. 

To see this observe that for all sufficiently small e > 0, R fe \ V is defined by 

£x 4 2 (^-l) 2 <e. 
i=i 

and has 2 fe connected components, since it retracts onto the set {0,l} fe . It now 
follows by Alexander duality that 

b k ^(V)=b (R k \V) = 2 k . 

Remark 1.2. Even though Theorem 11.11 is stated for semi-algebraic sets defined 
by a conjunction of weak inequalities, there is an easy reduction to this case for 
basic semi-algebraic sets defined by equalities and strict inequalities. The same 
reduction is also applicable to the algorithmic results described later in this paper, 
in particular Theorem 16.21 
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Note that the definition of cohomology groups of basic semi-algebraic sets defined 
by equalities and strict inequalities over an arbitrary real closed field R requires 
some care, and several possibilities exist. In this paper we follow the definition 
given in [9] which agrees with singular cohomology in case R = R. We refer the 
reader to [9] for further elaboration on this point. 

Wc have the following easy corollary of Theorem 11.11 

Corollary 1.3. Let R a real closed field and S C R fc be defined by 

f\ P = f\ P>0 f\ P < 

Pert Pev 2 Pev 3 

with deg(P) < 2 for each P <G Ui=o,i.2^' an ^ #Ui=o,i,2^ = s - 
Then, for all £ > 0, 

In the following proof, as well as later in the paper, we will extend the ground 
field R by infinitesimal elements. We denote by R(£) the real closed field of algebraic 
Puiseux series in £ with coefficients in R (see [5] for more details). The sign of a 
Puiseux series in R(£) agrees with the sign of the coefficient of the lowest degree 
term in £. This induces a unique order on R(£) which makes £ infinitesimal: £ is 
positive and smaller than any positive clement of R. When a £ R(£) is bounded 
from above and below by some elements of R, lim^(a) is the constant term of a, 
obtained by substituting for £ in a. Given a semi-algebraic set S in R fc , the 
extension of S to R', denoted Ext(S l , R'), is the semi- algebraic subset of R' defined 
by the same quantifier free formula that defines S. The set Ext(S', R') is well defined 
(i.e. it only depends on the set S and not on the quantifier free formula chosen to 
describe it) . This is an easy consequence of the transfer principle (see for instance 

Proof of Corollary \1.3[ Let < i C e < 1 be infinitesimals. We first replace 
the set S by the set S' C R{e) h defined by S' = Ext(S, R(e)) n B k (0, 1/e), where 
Bk(0,r) denotes the closed ball of radius r centered at the origin. It follows from 
Hardt's triviality theorem for semi-algebraic mappings [19] that bi(S) = bi(S') for 
all i > 0. We then replace S' by the set S" C R{e,6) k defined by, 

f\ P<0A-P<0 f\ -P + S<0 f\ -P-S <0 /\e 2 (Xl + ---+X%)-l<0. 

peVi Pev 2 Pev 3 

It follows from Hardt's triviality again that, bi(S') = bi(S") for all i > 0. Now 
apply Thcorcm ll.il □ 

Remark 1.4. Another point to note is that the bound in Theorem 11.11 depends 
crucially on the structure of the formula defining the semi-algebraic set S - namely 
that S is defined by a conjunction of polynomial inequalities of the form P < with 
deg(P) < 2. The polynomial bound on the highest Betti numbers no longer holds 
if S is defined by a formula involving quadratic inequalities but having a different 
structure - for instance, a disjunction instead of a conjunction of such inequalities. 
However, the polynomial bound continues to hold for certain small variations of 
the structure of the formula defining S. For instance, wc have the following slight 
generalization of Theorem 11.11 whose proof mimics that of Theorem 11.11 and which 
we omit. 
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Theorem 1.5. Let R a real closed field and ScR be defined by 

A ( V P « ^ 0), dcg(Py ) < 2, 1 < i < s, 1 < j < m. 

l<i<s 1 <J <rn 

Then, for £ > 0, 

b k - e (S)<(Sj\k°W. 

For fixed m and £, the above bound is still polynomial in s and k. Similarly, the 
main result of this paper, namely Algorithm |4] in Section [8j can be extended to this 
situation as well with complexity polynomial in s and k for fixed m and I. We will 
omit the details of this extension. 

Semi-algebraic sets defined by a system of quadratic inequalities have a special 
significance in the theory of computational complexity. Even though such sets might 
seem to be the next simplest class of semi-algebraic sets after sets defined by linear 
inequalities, from the point of view of computational complexity they represent 
a quantum leap. Whereas there exist (weakly) polynomial time algorithms for 
solving linear programming, solving quadratic feasibility problem is provably hard. 
For instance, it follows from an easy reduction from the problem of testing feasibility 
of a real quartic equation in many variables, that the problem of testing whether 
a system of quadratic inequalities is feasible is NPR-complete in the Blum-Shub- 
Smale model of computation (see [10])- Assuming the input polynomials to have 
integer coefficients, the same problem is NP-hard in the classical Turing machine 
model, since it is also not difficult to see that the Boolean satisfiability problem can 
be posed as the problem of deciding whether a certain semi-algebraic set defined 
by quadratic inequalities is empty or not (see Section [9|). Counting the number of 
connected components of such sets is even harder. In fact, we prove (see Theorem 
19. ip that for I = 0(log k), computing the £-th Betti number of a basic semi-algebraic 
set defined by quadratic inequalities in R fc is #P-hard. Note that PSPACE-hardness 
of the problem of counting the number of connected components for general semi- 
algebraic sets were known before [Ml [26] , and the proofs of these results extend 
easily to the quadratic case. In view of these hardness results, it is unlikely that 
there exist polynomial time algorithms for computing the Betti numbers (or even 
the first few Betti numbers) of such a set. In contrast to these hardness results, the 
polynomial bound on the top Betti numbers of sets defined by quadratic inequalities 
gives rise to the possibility that these might in fact be computable in polynomial 
time. 

In this paper we prove that for each fixed I > 0, the top £ Betti numbers of basic 
semi-algebraic sets defined by quadratic inequalities are computable in polynomial 
time. We will assume that the polynomials given as input to our algorithms have 
coefficients in some ordered domain D contained in a real closed field R. We will 
denote the algebraic closure of R by C. By complexity of our algorithms we will 
mean the number of arithmetic operations including comparisons in the ring D. 
When D = Z, we will also count the number of bit operations. 

The main result of this paper is the following. 
Main Result: We present an algorithm (Algorithm [4] in Section [8]) which given a 
set of s polynomials, V = {Pi, . . . , P s } C R[Ai, ...,X k ], with deg(P) < 2, 1 < i < 
s, computes bk-i{S), . . . , bk-e(S), where S is the set defined by Pi < 0, . . . , P s < 0. 
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The complexity of the algorithm is 

(i.i) e(:)^°" m 

If the coefficients of the polynomials in V are integers of bitsizes bounded by r, 
then the bitsizes of the integers appearing in the intermediate computations and 
the output arc bounded by r(sk) 2 ( ( ' )( . 

To our knowledge this is the first polynomial time algorithm for computing a 
non-trivial topological invariant of semi-algebraic sets in R fe defined by polynomial 
inequalities, where the number of constraints is not fixed and the polynomials are 
allowed to have degree greater than 1. The special case when all the polynomials 
are linear reduces to the well-studied problem of linear programming. In this case 
the set S is either a convex polyhedron or empty, and (weakly) polynomial time 
algorithms are known to decide emptiness of such a set. In another direction, 
Barvinok [2] designed a polynomial time algorithm for deciding feasibility of systems 
of quadratic inequalities, but under the condition that the number of inequalities 
is bounded by a constant (see also [17] for an interesting generalization as well as 
a constructive version of this result). 

2. Brief Outline 

Given any compact semi-algebraic set 5 1 , we will denote by bi(S) the rank of 
H l (S,<Q) (the i-th simplicial cohomology group of S with coefficients in Q). We 
denote by S k C R fe+1 the unit sphere centered at the origin. We first consider the 
case of semi-algebraic subsets of the unit sphere, S fc C R fc+1 , defined by homoge- 
neous quadratic inequalities. We then show how to reduce the general problem to 
this special case. 

Let S C S fc be the set defined on S fc by s inequalities, Pi < 0,...,P S < 0, 
where Pi, . . . , P s £ R[Vo, ■ ■ ■ , -Xfc] are homogeneous quadratic polynomials. For 
each i, 1 < % < s, let 5, C S fc denote the set defined on S k by P t < 0. Then, 
S = nf—jSj. There are two main ingredients in the polynomial time algorithm for 
computing the top Betti numbers of S. 

The first main idea is to consider S as the intersection of the various S^s and 
to utilize the double complex arising from the generalized Mayer- Vietoris exact 
sequence (see Section [3]). It follows from the exactness of the generalized Mayer- 
Vietoris sequence (see Proposition [1] below) , that the top dimensional homology 
groups of S are isomorphic to those of the total complex associated to a suitable 
truncation of the Mayer- Vietoris double complex. However, computing even the 
truncation of the Mayer- Vietoris double complex, starting from a triangulation 
of S would entail a doubly exponential complexity. However, we utilize the fact 
that terms appearing in the truncated complex depend on the unions of the Si's 
taken at most £ + 2 at a time. There are at most if) such sets. Moreover, 

for semi-algebraic sets defined by the disjunction of a small number of quadratic 
inequalities, we are able to compute in polynomial (in k) time a complex, whose 
homology groups arc isomorphic to those of the given sets. The construction of these 
complexes in polynomial time is the second important ingredient in our algorithm 
and is described in detail in Section [7] ( Algorithm [2]) . These complexes along with 
the homomorphisms between them define another double complex whose associated 
spectral sequence (corresponding to the column-wise filtration) is isomorphic from 
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the E2 term onwards to the corresponding one of the (truncated) Mayer- Vietoris 
double complex (see Theorem 16.21 below) . Since, we know that the latter converges 
to the homology groups of S, the top Betti numbers of S are equal to the ranks 
of the homology groups of the associated total complex of the double complex we 
computed. These can then be computed using well known efficient algorithms from 
linear algebra. 

The rest of the paper is organized as follows. In Section [3] we recall certain 
basic facts from algebraic topology including the notions of complexes, and double 
complexes of vector spaces, spectral sequences and triangulations of semi-algebraic 
sets. We do not prove any results since all of them are quite classical and we refer 
the reader to appropriate references p~2l [211 E] for the proofs. In Section 21 we 
recall some basic algorithms in semi- algebraic geometry that we will need later. 
We state the inputs, outputs and complexities of these algorithms. For classical 
results in the field of real algebraic geometry, as well as the details of certain 
algorithms that we use, we refer the reader to [9] for convenience. In Section [5] we 
describe certain topological properties of semi-algebraic sets defined by quadratic 
inequalities which are crucial for our algorithm. Most of the results in this section 
are due to Agrachev [1] . In Section [6] we prove the main mathematical results 
necessary for our algorithm. In Section [7] we describe our algorithm for computing 
the top Betti numbers of semi-algebraic sets defined by quadratic forms. We treat 
the general case in Section [8] Finally, in Section [9l we show the computational 
hardness of computing the first few Betti numbers of a given semi-algebraic set 
defined by quadratic inequalities, by proving that the problem is #P-hard. 

3. Topological Preliminaries 

In this section we recall some basic facts from algebraic topology, related to 
double complexes, and spectral sequences associated to double complexes. We also 
fix our notations for these objects. All facts that we need are well known, and we 
merely give a brief overview, referring the reader to [T2l [2Tj for detailed proofs. 

3.1. Complex of Vector Spaces. A sequence {C p }, p € Z, of Q-vector spaces 
together with a sequence {d p } of homomorphisms d p : C p — -> C p+1 for which 
QP+ 1 o d p = for all p is called a complex. 

The cohomology groups, H P (C',Q) are defined by, 

H p (C',Q) = Z p {C')/B p {C), 

where B P {C) = Im(d p ~ r ), and Z P {C) = Kcr(5 p ). 

The cohomology groups, H*(C',F), are all Q-vector spaces (finite dimensional 
if the vector spaces C p 's are themselves finite dimensional). We will henceforth 
omit reference to the field of coefficients Q which is fixed throughout the rest of the 
paper. 

Given a complex C* , we denote by C. the dual complex, 

< <^p-i < <-'p < 

where C p = Hom(C p ,Q) is the vector space dual to C p and 

d p : Hom(C p ,Q) -> Hom(C p_1 ,Q) 

is the homomorphism dual to d 13 " 1 . 
Moreover, H*{C„Q) = H*(C',Q). 
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Given two complexes, C* = (C p ,d p ) and D* = (D p ,d p ), a homomorphism of 
complexes, <j> : C* — > is a sequence of homomorphisms ^ : C p — > Z) p for which 
d p o(j) p = (t> p+1 o 9 p for all p. 

In other words, the following diagram is commutative. 



C p 
D p 



C p+i _ 



A homomorphism of complexes, <j) : C* — > D", induces homomorphisms, </>* : 
H*(C) — > H*(D'). The homomorphism cj> is called a quasi- isomorphism if the 
homomorphisms 0* arc isomorphisms. 



3.2. Double Complexes. A double complex is a bi-graded vector space, 

C'' = C p ' q , 

with co-boundary operators d : C p ^ -> C^ 1 and <5 : C p > q -> C p+1 '« and such that 
e?(5 + <5c? = 0. We say that C*'* is a first quadrant double complex, if it satisfies the 
condition that C p ' q = if either p < or q < 0. Double complexes lying in other 
quadrants are defined in an analogous manner. 



d 

C 0,2 c 



d 

a 

d 

& 



o,i S M S 



c 1 



c 1 



a* 



,1 s 



We call the complex, Tot*(C"'*), defined by 

Tot n (C"'*) = C p ' q , 

p+q=n 

ith differential 



D™ = d + S : Tot^C*'*) -> Tor+^C"'*) 
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to be the associated total complex of C*. 





3.3. Spectral Sequences. A spectral sequence is a sequence of bigraded complexes 
(E r , d r : Epi -> EP +r ' q - r+1 ) (see Figure|T]) such that the complex E r+1 is obtained 
from E r by taking its homology with respect to d r (that is E r +i — Hd r (E r )). 

There are two spectral sequences, 'E^' q ,"E^' q , (corresponding to taking row- 
wise or column- wise nitrations respectively) associated with a double complex C*'*, 
which will be important for us. Both of these converge to if* (Tot* (C*'*)). This 
means that the homomorphisms d r are eventually zero, and hence the spectral 
sequences stabilize, and 

>E™= "^=^(Tof(C-)), 



p+q=i 



p+q=i 



for each i > 0. 

The first terms of these are 'E x = H S {C '),'E 2 = H d H s (C''), and "E 1 = 
H d (C->'),"E 2 =H s H d (C'> a ). 

Given two (first quadrant) double complexes, C*'* and C*'*, a homomorphism 
of double complexes, 



is a collection of homomorphisms, <j> p,q : C pq — 
diagrams commute. 



C p ' q , such that the following 



C p ' q 



i,p,q 
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(JP-q d y (jp,q+l 



(jp q d ; (jp.q+i 
A homomorphism of double complexes, 

<t> : C' u — ► C'' m , 

induces homomorphisms, <f) s : E s — > E s between the associated spectral sequences 
(corresponding either to the row- wise or column- wise nitrations). For the precise 
definition of homomorphisms of spectral sequences, see [H]. We will use the fol- 
lowing useful fact (see [2T], page 66, Theorem 3.4) several times in the rest of the 
paper. 

Theorem 3.1. If <p s is an isomorphism for some s > 1, then E^ q and E^ q 
are isomorphic for all r > s. In other words, the induced homomorphism, <j) : 
Tot*(C* : *) — > Tot*(C* : *) is a quasi-isomorphism. 

3.4. Triangulation of semi-algebraic sets. A triangulation of a compact semi- 
algebraic set S is a simplicial complex A together with a semi-algebraic homeomor- 
phism from | A| to S. Given such a triangulation we will often identify the simpliccs 
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in A with their images in S under the given homeomorphism, and will refer to the 
triangulation by A. 

Given a triangulation A, the cohomology groups H l (S) are isomorphic to the 
simplicial cohomology groups H l (A) of the simplicial complex A and are in fact 
independent of the triangulation A (this fact is classical over R; see for instance [9] 
for a self-contained proof in the category of semi-algebraic sets). 

We call a triangulation hi : |Ai| — > 5* of a semi-algebraic set S, to be a refinement 
of a triangulation /12 : IA2I — > S if for every simplex o\ € Ai, there exists a simplex 
(T 2 E A 2 such that h\(oi) C /^(o^)- 

Let Si C 5 2 be two compact semi-algebraic subsets of R fe . We say that a semi- 
algebraic triangulation h : |A| — > S2 of S2, respects Si if for every simplex a G A, 
h(a) fl Si — h(a) or 0. In this case, h~ 1 (Si) is identified with a sub-complex of A 
and h\h-i(Si) '■ h~ 1 (Si) — > Si is a semi-algebraic triangulation of Si. We will refer 
to this sub-complex by A^. 

We will need the following theorem which can be deduced from Section 9.2 in 
PU (see also 0). 

Theorem 3.2. Let Si C S2 C R fc 6e closed and bounded semi-algebraic sets, and 
let hi : Ai Si, i = 1,2 be semi- algebraic triangulations of Si,S2- Then, there 
exists a semi- algebraic triangulation h : A — ► 5*2 of S2, such that A respects Si, A 
is a refinement of A2, and A\s 1 is a refinement of Ai. 

Moreover, there exists an algorithm which computes such a triangulation whose 
complexity is bounded by (sd) 2 ( , where s is the number of polynomials used in 
the definition of Si and S2, and d a bound on their degrees. 

3.5. Mayer- Vietoris Double Complex. Let Si,...,S s C R fe be closed and 
bounded semi-algebraic sets, and let S = n* =1 Sj. Choose a sufficiently fine trian- 
gulation of Uf =1 5<, such that all intersections of the form, Si (l - ••r\Si t , correspond 
to subcomplexes of the simplicial complex of the triangulation. Note that the exis- 
tence of such a triangulation (in fact, a semi- algebraic triangulation) is well known 
(see 0). However, the best algorithm for computing such a triangulation has com- 
plexity which is doubly exponential (in k) and produces a doubly exponential sized 
triangulation, and hence is not suitable for our purpose. 

The first main ingredient for our polynomial time algorithm is the double com- 
plex associated to the generalized Mayer- Vietoris sequence, which we describe first. 

For each Si, let Ai be the subcomplex corresponding to it, and let A = Ai (1 ■ ■ ■ D 
A s . For any p > 0, let A a °----^ denote the union, A ao U- • -UA ap . Let C t {A) denote 
the Q- vector space of i-chains of A, and C,(A) the corresponding chain complex. 

Proposition 1. The following sequence is exact. 

O^C,(4)A0C,(r)^ C.(A ao ' ai ) ••• 

-L> c .(A ao >"> a *) -$-> C.(A a °<"' a >+ 1 ) ••• , 

OD<-<a p Q <---<a p+ i 

where i is induced by inclusion and the connecting homomorphisms 5 arc defined 
as follows: 

for C € ® C[0 <•■■<a p C.(^4 Qo '•••' Q!, ), (Sc) a0t ,„ t a p+1 = J2 < i < p+ i(~lYCao,...,d i ,...,a p+1 - 

Proof. See [5]. □ 
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We have a corresponding (fourth quadrant) double complex A/"*'*, with 







N v ' q 




■ ,a p \ 
























— > 




6 ; 


©QQ<Ql Cfc {A ' 1 ) 


S 












— > 


® ao C k ^(A ao ) 


6 


$QO<ai Cfc — 1 (-<4 ° ' 1 ) 


5 




l» 








— > 


© Q0 C* fe _ 2 (A a °) 






5 




1" 




l» 




— > 


© Q0 C fe _ 3 (A Q0 ) 




® ao<ai C k - 3 (A a °^) 


5 } 




1" 











l« 

\±Jao<ai <Q2 
ao<a 1 <a 2 { - y k — 2\-fi J 

l» 

7tt <tti<a 2 t - / fc-3l /1 " J 

1" 



Standard facts about the Mayer- Vietoris spectral sequence then yields the fol- 
lowing well known proposition (see [S]). 

Proposition 2. For < i < k, 

Hi(A) £ JP(Tbr (Af •*)). 

Note that since A/"*'* is a fourth-quadrant double complex, the complex Tot" (TV*'* ) 
is defined by, 

Tot^/V*'*) = /V p ' 9 . 

p+k—q—i 

Moreover, if we denote by M*'* the truncated complex defined by, 

Mf q = N p ' q , Q<p<l+l, k-i-l<q<k, 
= 0, otherwise, 

then it is clear that, 

(3.1) Hi{A) ^H^Totf (M''*)), k~£<i<k. 

As noted previously, we cannot hope to compute even the truncated complex 
A/7'* since we do not know how to compute triangulations efficiently. We overcome 
this problem by computing another double complex 2?*'*, such that there exists a 
homomorphism of double complexes, 

V> : ^ A/7'\ 

which induces an isomorphism between the E\ terms of the spectral sequences asso- 
ciated to the double complexes T>*' % and A/J'*. This implies, by virtue of Theorem 
13.11 that the homology groups of the associated total complexes are isomorphic, 
that is, 

H*(Tof(Af''')) S iT (Tot* (X>*<*)). 

The construction of the double complex T>\'' is described in the Section [6] 
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4. Algorithmic Preliminaries 

In this section, we recall some basic algorithms from semi-algebraic geometry 
that we will need later in the paper. Our main reference is [5]. Here we just 
describe the inputs, outputs and the complexities of these algorithms. 

4.1. Sign Conditions. Let R be a real closed field and C its algebraic closure. 
A sign condition is an element of {0, 1, —1}. We define 

sign(.T) =0 if and only if x — 0, 
\ sign(cc) = 1 if and only if x > 0, 
sign(.T) = — 1 if and only if x < 0. 

Let Q C R[Ai, . . . , Xk]. A sign condition on Q is an element of {0,1, -l} 9 . We 
say that Q realizes the sign condition a at x € R fe if 

/\ Bign(Q(aO) = <t(Q). 
QeQ 

The realization of the sign condition a is 

11(a) = {x e R k | f\ sign(Q(x)) = a(Q)}. 
QeQ 

The sign condition a is realizable if 1Z(a) is non-empty. We denote by sign(Q) the 
set of realizable sign conditions of Q. 

4.2. Representations of points. Since our algorithms will have to deal with 
points whose co-ordinates are not rational numbers, we first describe the particular 
representations of such points that we will use. In the following D denotes an 
ordered domain contained in the real closed field R. 

Let P E R[X] and a £ {0, 1, -ljEM^ a s i gn condition on the set Der(P) of 
derivatives of P. The sign condition a is a Thorn encoding of x £ R if a(P) = 
and a is the sign condition taken by the set Der(P) at x. A real root x of P is 
completely characterized by the signs of the derivatives of P at x and we say that 
x is specified by a. Given a Thorn encoding er, we denote by x(a) the root of P in 
R specified by a. 

The representations of points are as follows. A fc-univariate representation is a 
k + 2-tuple of polynomials of D[T], 

(f(T),g (T), gi (T),...,g k (T)), 

such that / and go are coprime. Note that go(t) ^ if t £ C is a root of f(T). The 
points associated to a univariate representation are the points 

( ffiW 9kW \ rfe 

\go(t)""'go(t)) 

where t € C is a root of f(T). 

A real ^-univariate representation a pair u, <r where u is a fc-univariate represen- 
tation and a is the Thorn encoding of a root of /, t a G R. The point associated to 
the real univariate representation is the point 

9i(t<r) 9k(t a ) \ g Rfc 
9o(U) ' " ' ' 9o(ta) J 
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4.3. Complexity of Linear Algebra Subroutines. In our algorithms, we will 
also have to compute eigen-vectors of, as well as perform Gram-Schmidt orthogo- 
nalizations on, real symmetric matrices whose entries are described by univariate 
representations. Suppose that M is a matrix of size k x k, given as a real k 2 - 
univariatc representation. Recall that this means that M is given as a pair, cr, u, 
where cr is a Thorn encoding and u = (/(T), go(T), gi i i(T), gk,k{T)) is a k 2 - 
univariatc representations. 

If the degrees of the polynomials /, <?o> <?i,i) • • • j <?fc,fc are bounded by D, then 
using the most naive algorithms, all the above tasks can be performed with com- 
plexity bounded by (kD)°^ . In order to see this notice that the complexity of 
the algorithms of linear algebra with matrices of size k x k is bounded by k ^ 1 ' . 
However, since the computations takes place in the ring D[T]/(/), each arithmetic 
operation in this ring costs operations in the ring D. 

4.4. Complexity of Computing Linear Arrangements. We will also need to 
compute descriptions of cell complexes, whose cells are the chambers of different di- 
mensions in an arrangement of £ hypcrplanes in R fe . Notice that the number of cells 
in such an arrangement can be bounded naively by 3 e . Using standard algorithms, 
descriptions of the cells in such a complex, and their incidence relationships, can 
be computed with complexity 

Moreover, in case the arrangement is parametrized, that is the hypcrplanes arc 
described by linear equations, whose coeffcients are described by parametrized uni- 
variate representations, with at most I parameters, and whose degrees are bounded 
by D, one can compute a family of polynomials in the parameters, such that for 
all values of the parameters satisfying a fixed sign condition on this family, the 
combinatorial structure of the cell complex remain constant. Moreover, the com- 
plexity of this procedure, as well as the number and degrees of this new family of 
polynomials, are all bounded by {kD)°^K 

5. Topology of sets defined by quadratic constraints 

The results of this section were proved by Agrachev pQ in the context of non- 
degenerate (see [I] for the precise definition of non-dcgcncracy) quadratic maps. 
However, for the purposes of this paper the assumption of non-degeneracy is not 
required as shown below. 

Let Pi, . . . ,P S be homogeneous quadratic polynomials in R[Xo, . . . , Xk\. 

We denote by P = (Pi , . . . , P s ) : R fe+1 -> R s , the map defined by the polynomials 

P,..., P. 
Let 

T= |J {x G S fc I Pi(x) < 0}. 

l<i<s 

Let 

n = {uj e R s | M = < 0, 1 < i < s}. 
For w 6 wc denote by uiP the quadratic form defined by 

s 

up = Y,"iPi- 

i=l 

Let B C O x S fe be the set defined by, 

B = {(uj, x) I uj G Q,, x G S fc and ojP(x) > 0}. 
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We denote by 4>\ : B — > il and </>2 : B — > S k the two projection maps. 



B 




n s 



With the notation developed above, 

Proposition 3. The map (f>2 gives a homotopy equivalence between B and 4>2{B) = 
T. 

Proof. We first prove that 02 (-B) = 7 1 . If ir G T, then there exists some £, 1 < i < s, 
such that Pi(x) < 0. Then for w = (— Sij, . . . , — 8 St i) (where 5ij = 1 if i = j, and 
otherwise), we see that (w,x) £ £?. Conversely, if x e 412(B), then there exists 
cj = (ct>i, . . . ,u s ) € fi such that, J^i=i ^P^) ^ 0- Since, < 0, 1 < i < s, and 
not all Wj = 0, this implies that Pi(x) < for some i, 1 < i < s. This shows that 
x £ T. 

For x G <j) 2 (B), the fibre 

^2 1 (a:) = {(w, x) I u G f2 such that wP(x) > 0}, 

is a non-empty subset of Q defined by a single linear inequality. Thus, each non- 
empty fiber is an intersection of a convex cone with S'\ Thus, all such fibres can 
be retracted to their respective centers of mass continuously, proving the first half 
of the proposition. □ 

For any quadratic form Q, we will denote by index(Q), the number of negative 
eigenvalues of the symmetric matrix of the corresponding bilinear form, that is of 
the matrix M such that, Q(x) = (Mx,x) for all x G K k+1 (here (•, •) denotes the 
usual inner product). We will also denote by Aj(Q),0 < i < k, the eigenvalues of 
Q, in non-decreasing order, i.e. 

A (Q) < Ai(Q) < ••■< A fe (Q). 

Given a quadratic map P = (Pi, . . . , P s ) : R k+1 -> R s , and < j < k, we denote 

by 

Qj ■ = {lo e n I Xj(uP) > 0}. 
It is clear that the fij's induce a filtration of the space f2, i.e., fio C fii C • • • C 

Agrachev [I] showed that the Leray spectral sequence of the map <j>\ (converging 
to the cohomology H*(B) = H*(T)), has as its E 2 terms, 

(5.1) £f = fP(n fc _ g ,n fc _ g _i). 

(We refer the reader to [Hj for the precise definition of the Leray spectral sequence 
of a map). 

Equation 15.11 follows directly from the following lemma. 

Lemma 5.1. The fibre of the map 0i over a point u] G Clj \ has the homotopy 

type of a sphere of dimension k — j. 
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Proof. First notice that for co £ Qy \ Clj—i, \o(u)P), . . . , Xj-i(u>P) < 0. Moreover, 
letting Yq(ujP), . . . ,Yk(uiP) be the co-ordinates with respect to an orthonormal 
basis consisting of the eigenvectors of wP, we have that is the subset of S 

defined by, 

k 

^2\ 1 {ujP)Y 1 {ujP) 2 > 0, 

i=0 

Y,Yi^Pf = l. 

Since, Ai(wP) < 0, < i < j, it follows that for co 6 flj \ the fiber 

(j)^ 1 ^) is homotopy equivalent to the (fc — j)-dimcnsional sphere defined by setting 
Y q (ujP) = ■■■ = Yj-i(ujP) = on the sphere defined by £)* =0 Y l {toP) 2 = 1. □ 

6. Computing the cohomology groups of a basic semi-algebraic set 
defined by homogeneous quadratic inequalities 

In this section, we will show how to effectively compute the spectral sequence 
described in the previous section. 

Let V = (Pi, . . . , P s ) C R[^07 ■ ■ ■ j Xk] be a s-tuple of quadratic forms. For any 
subset Q C P, we denote by Tg C S fc , the semi-algebraic set, 

T Q = [}{xe S k I P(x) < 0}, 
-Pee 

and let 

S = p| {x e S k | P(x) < 0}. 

We denote by C(H(Tq)) the co-chain complex of a cellular subdivision, H(Tq) of 
Tg, which is to be chosen sufficiently fine (to be specified later). 

We first describe for each subset Q CP with #Q = I < fc, a complex, and 
natural homomorphisms, 

i> Q :C'(H(T Q ))^M' Q , 

which induce isomorphisms, 

r Q :H*(C'(H(T Q )))^H*(M' Q ). 

Moreover, for B C A C P with j^A = < k, we construct a homomorphism 

of complexes, 

4>A,B ■ M' A -> M%, 
such that the following diagram commutes, 



(6.1) 



H*(M' A ) 



U A,B 



H*(M' B ) 



H*(C'(H(T A ))) — H*(C'(H(T B ))) 
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where (\>\ B and r* are the induced homomorphisms of 4>a..B and the restriction 
homomorphism r respectively. 

Now, consider a fixed subset Q c P, which without loss of generality we take to 
be {P u ...,P e }. Let 

Pq = (P,...,P,) :R fe+1 ^R £ 
denote the corresponding quadratic map. 

As in the previous section, let R® = R £ , and 

fig = {uj £ R e | |cj| = 1, uji < 0, 1 < i < £}. 

Let Bq C fig x S fc be the set defined by, 

B Q = {(^i) I W G Qq,x S S fc andwPgfz) > 0}, 

and we denote by <f>i q : Pg — > fig and </>2,g : Bq — > S fc the two projection maps. 
For each subset Q' C <2 we have a natural inclusion fig' fig. 

6.1. Index Invariant TViangulations. We now define a certain special kind of 
semi-algebraic triangulation of fig that will play an important role later. 

Definition 6.1. (Index Invariant Triangulation) An index invariant triangulation 
of fig consists of: 

(A) A semi-algebraic triangulation, 

: — > fig 

of fig, which is compatible with the subsets fig/ for every Q' C Q, and 
such that for any simplex a of Ag, index(wPg), as well as the multiplicities 
of the eigenvalues of ujPq, stay invariant as u> varies over /iq(<t), and 

(B) for every simplex a of Ag, with indexfw Pg) = i, a continuous choice of 
an orthonormal basis of R fe+1 , {eo(cr, w), . . . , efc(cr, w)}, for a; £ ^g(c), such 
that {ej(cr, w), . . . , e-k{a, uj)} span the linear subspace of R fc+1 on which the 
quadratic form loPq is positive semi-definite. 

Wc describe later how to compute index invariant triangulations (see Algorithm^ 
below). It will follow from the complexity analysis of Algorithm[T]that the size of the 
complex Ag as well as the degrees of the polynomials occuring in the parametrized 
univariate representations defining the various bases, {eo(cr, w), . . . , efe(<r, to)}, are 
all bounded by fc 2 ° w . 

For the rest of this section we fix an index invariant triangulation /ig : Ag — > fig, 
satisfying the complexity estimates stated above. 

The following proposition relates, for any simplex a £ Ag, the homotopy type 
of </> 1 ~g(/ig(c)) to that of a single fiber. 

Proposition 4. For any simplex a £ Ag and uj £ /ig(cr), 4 ) iq(1 1 Q.( ~)) is homotopy 
equivalent to g(w), and both these spaces have the homotopy type of the sphere 

g/c— index(u;P) 

Proof. Let i = index(wP). Since indcxfwP) is invariant as u> varies over /ig(cr), 
the quadratic forms uiP has exactly i negative eigen- values for each uj £ ft,g(cr). Let 
M(a, uj) C R fe+1 be the orthogonal complement to the linear span of the correspond- 
ing eigen- vectors, and let B(a, uj) = M(a, w)flS fe . Clearly, M(cr, uj) and P(cr, uj) vary 
continuously with u>, and <^>j~g(w) can be retracted to the set {uj} x B(a, uj). Finally, 
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since /iq(ct) is contractible to u>, its clear that </>! q(/iq(ct)) retracts to {uj} x B(a, u) 
and the latter has the homotopy type of s fc_indcx(wP) by Lemma IO □ 

6.2. Definition of the cell complex JC(Bq). Our next goal is to construct a cell 
complex homotopy equivalent to Bq obtained by glueing together certain regular 
cell complexes, /C(er), where a £ Aq. 




Figure 2. The complex Aq. 



C,nC„,,nC„,, 




Figure 3. The corresponding complex C(Ag). 

Let 1 > eo > £1 > ■ • ■ > e s > be infinitesimals. For r\ S Aq, we denote by 
C,, the subset of fj defined by, 

C v = {x s fj | and dist(x, 6*) > £dim(0) f° r au $ ~< f 7 }- 

Now, let cr -< rj be two simplices of Aq. We denote by C a ^ the subset of 77 
defined by, 

C a . n = {x E fj I dist(x, a) < £dim(<r)j an d dist(x, 9) > £dim(0) for all 9 -< a}. 
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Note that, 

|a q |= (J au |J C a „. 

creAg tr,r;£Ag,cr-;?) 

Also, observe that the various C v 's and (7 CTlJ7 's are all homeomorphic to balls, and 
moreover all non-empty intersections between them also have the same property. 
Thus, the union of the C^'s and C a ^s together with the non-empty intersections 
between them form a regular cell complex, C(Aq), whose underlying topological 
space is |Ag| (see Figures [2] and [3|) . 

We now associate to each C CT (respectively, C a:V ) a regular cell complex, K(a), 
(respectively, K,{a,rf)) homotopy equivalent to (^^(^(Co-)) (respectively, 

For each a G Ag, and lo G /iq(c), let 

A CT M = • • • = A£(w) < Xl +1 H = • • • = Af» < • • • < = • • • = 

Af p M < o = ••• = Af p+1 M <... = ... = a^h 

denote the eigenvalues of toP. Here, index(tjP) = i p + 1. Also, since the mul- 
tiplicities of the eigenvalues do not change as lo varies over /ig(er), the block 
structure, [0, . . . , io], [io + 1, . . . , ii], ...,[-,..., fc] also does not change as lo varies 
over /iq(<t). For < j < p, let M^{a,Lo) denote the subspace of R fc+1 orthogo- 
nal to the subspace spanned by the eigenvectors corresponding to the eigenvalues 

KM = ••• = KM < K 0+ M = ••• = KM < ••• < K p - X+ M = ••• = 

Af. (u>), and let M(cr, w) = M p (a, lo). Since the eigenvalues vary continuously and 
their multiplicities do not change as lo varies over /iq(<t), the flag of subspaces 
M°((T,lo) D ■ ■ ■ D M p (<j,uj) also varies continuously over /ig(cr). 

For each a G Ag, and lo G /iq(ct), let {eo(cr, w), . . . , efc(cr, w)}, be the continuously 
varying orthonormal basis of R fc+1 computed previously. 

We extend the the orthonormal basis, {eo(cr, lo), . . . , efc(cr, lo)}, continuously to 
each C a _ n for rj with a -< 77, satisfying the condition that 

M(a,uo) C span(e 4 (cr,w), . . . ,e fc (cr,o))). 

This extension can be done in a consistent manner because the first i eigenvalues, 
Ao(wP), . . . , Aj_x(wP) of uP stay negative, and Ai_i(wP) < Xi(wP) for w in any 
infinitesimal neighborhood of hQ{C a ). Thus, the linear subspace of R fc orthogo- 
nal to the eigenspaces corresponding to the eigenvalues, Xq(loP), . . . , Aj_i(c<;P) is 
well defined, and varies continuously with lo in any infinitesimal neighborhood of 
h Q {C a ). 

More precisely, for any point z' = /ig 1 (cj') G C CT)?) , let a; = }iq(lo), be the unique 
point in C a closest to z' . For j = i, . . . , k, let e'j(cr, lo') be the orthogonal projection 
of ej(<T,ui) onto the subspace M(a,Lo'), and let ei(a,u>'), . . . , ek(cr, lo') be obtained 
from e^(a, u/), . . . , e' k (cr, lo') by Gram-Schmidt orthogonalization. Note that each 
e.j{<J, lo') obtained this way is infinitesimally close to ej(a,Lo). Using the same 
procedure starting with the vectors eo(<j,u>), . . . ,ei_i(er, lo) and the linear subspace 
M(a,oj') ± , we compute an orthonormal set ofvectors eo(a,Lo'), . . . , ej_i (a, lo') span- 
ning M{a, lo')- 1 -. It is clear from construction, that the vectors eo(<r,Lo'), . . . , efc(er, lo') 
form an orthonormal basis for R fe+1 . Moreover, they vary continuously with lo', 
extending the basis eo(er, lo), . . . , efc(er, lo), and finally ej(cr, lo'), . . . , efc(er, lo') span 
M(a,Lo'). 
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The orthonormal basis 

{e (a,(j), . . .,e k (a,uj)}, 

determines a complete flag of subspaces, J- (a, to), consisting of 

L°(<r,u) = 0, 
L 1 (cr,o;) = span(efc(er, w)), 
L 2 (a,uj) = span(e fc (cr,ai),efc_i(cr,w)), 

L k+1 (a 7 oj) = R fe+1 . 

For < j < k, let c+(cr,w) (respectively, Cj (a,uj)) denote the (k — j)-dimensional 
cell consisting of the intersection of the L fe ~ • 7+1 (c, u;) with the unit hemisphere in 
R fc+1 defined by {x € S k | (x, ej(a, u)) > 0} (respectively, {x G S k | (x,ej(a,(j)) < 
0}). 

The regular cell complex K,{a) (as well as K.(a,r])) is defined as follows. 

The cells of IC(a) are {(x,uj) | x S c^{ct,uj),u) € /ig(c)}, where index(wP) < j < 
k, and c £ C(Ag) is either C a itself, or a cell contained in the boundary of C a . 

Similarly, the cells of IC(a,r/) are {(x,w) | x € c^(cr,uj),u> € ^q(c)}, where 
index(o;P) < j < k, c S C(Ag) is either C ari itself, or a cell contained in the 
boundary of C a _ n . 

Our next step is to obtain cellular subdivisions of each non-empty intersection 
amongst the spaces associated to the complexes constructed above, and thus obtain 
a regular cell complex, IC(Bq), whose associated space, \K,(Bq)\, will be shown to 
be homotopy equivalent to Bq (Proposition [6] below). 

First notice that /C(cr', rj')\ (respectively, |/C(<r)|) has a non-empty intersection 
with |/C(cr, rf)\ only if C a i^i (respectively, C a >) intersects C aiT] . 

Let C be some non-empty intersection amongst the C^s and C<j ir; 's, that is C is 
a cell of C(Aq). Then, C C r\ for a unique simplex rj 6 Aq, and 

c = Ceri,i) n • • • n c aptri , 

with a\ < (j-2 -<■■■< cr p -< r] and p < #Q + 1. 
Consider u) £ ^-q(C). We have p different flags, 

J"((Ti,w), . . .,T{a p ,u), 

giving rise to p independent regular cell decompositions of B{lo, rj) = M(j], to) n S fc . 

There is a unique smallest regular cell complex, K.'(C,oj), that refines all these 
cell decompositions. The cells of this cell decomposition consists of the following. 
Let L C M(j],uj) be any linear subspace of dimension m, < m < k + 1, which is 
an intersection, of linear subspaces L\, . . . , L p , where Z, £ J-(ai,uj), 1 < i < p. The 
elements of the flags, T{p\, oj), . . . , T{a p , uj) of dimensions m + 1, partition L into 
polyhedral cones of various dimensions. The union of the sets of intersections of 
these cones with S fc , over all such subspaces L C M(rj,uj), are the cells of /C'(C, uj). 
Figure 2] illustrates the refinement described above in case of two flags in R 3 . 

We now triangulate /iq(C), using the algorithm implicit in Theorem l3.2l (Trian- 
gulation), such that the combinatorial type of the arrangement of flags, 

J r (a ll uj) 1 . . .,T{a p ,uj) 
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Figure 4. The cell complex K'(C,u). 

and hence the cell decomposition K'(C,U)), stays invariant over the image, hc(0), 
of each simplex, 6, of this triangulation. More precisely, we first compute a family 
of polynomials, Ac C R[Zi, . . . , Zg\ whose signs at to determine the combinatorial 
type of the corresponding arrangement of flags. It is easy to verify (see Section 
I4.4p . given the complexity bounds on the parametrized univariate representations 
defining the orthonormal bases, {eo(c, u>), . . . , efc(tr, u>)}, u> £ ^-q(o'), stated above, 
that the number and degrees of the polynomials in the family Ac is bounded by 
k 2 { ' . We then use the algorithm implicit in Theorem 13.21 (Triangulation) , with 
Ac as input, to obtain the required triangulation. 
The closures of the sets 

{(w,x) | xeceK'(C,u),ueha(ho(e))} 

constitute a regular cell complex, tC{C), which is compatible with the regular cell 
complexes fC(ai), . . . , IC(a p ). 

The following proposition gives an upper bound on the size of the complex K,(C). 
We use the notation introduced in the previous paragraph. 

Proposition 5. For each w £ Hq{C), the number of cells in fC'(C,u>) is bounded by 
k°( e \ Moreover, the number of cells in the complex /C(C) is bounded by k 2 ° H) . 

Proof. The first part of the proposition follows from the fact that there are at most 
£#Q+i _ choices for the linear space L and the number of (m — 1) dimensional 
cells contained in L is bounded by 2 l (which is an upper bound on the number of 
full dimensional cells in an arrangement of at most t hypcrplancs). The second part 
is a consequence of the complexity estimate in Theorem l3. 21 (Triangulation) and the 
bounds on number and degrees of polynomials in the family Ac stated above. □ 
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We denote by K(Bq), the union of all the complexes /C(C) constructed above, 
noting that by construction, K,(Bq) is a regular cell complex. 

Proposition 6. |/C(Bq)| is homotopy equivalent to Bq. 

Proof. We first show that Bq is homotopy equivalent to a subset Bq C Bq as 
follows. For each simplex a G Aq of the largest dimension £, we use the retraction 
used in the proof of Proposition^ to retract 1 " 1 (/iQ(relint(C (J ))) to the set {(w,x) | 
uj G relint(C (J ), x G B(cj,cr)}. In this way we obtain a semi-algebraic set, X' t , 
which is a deformation retract of Ext(£?g, R(eo, • • • > £i-i))- Let Xg = lim e ^_ 1 X' t 
Notice that in the definition of X' e , if we replace eg-\ by a variable t and denote the 
corresponding set by X' t t , then for all < i < t' , X' t t C X' e t , and each X^t is closed 
and bounded. It then follows (see Lemma 16.17 in [9]) that Ext(X^, R(eo, ■ ■ ■ , £i-i)) 
has the same homotopy type as XL and hence Xg has the same homotopy type as 
Ext(B g ,R(e , ■ ■ ■ ,££-2>)- 

Now repeat the process using the (£ — l)-dimensional simplices and so on, to 
finally obtain Xq = B'q, which by construction has the same homotopy type as Bq. 
Finally, (again using Lemma 16.17 in [S]) we also have that Xq = lim eo |/C(_Bg)| 
and Ext(Xo, R(eoj • ■ ■ , e i-i)) nas the same homotopy type as \K{Bq)\. □ 

We also have, 

Proposition 7. The number of cells in the cell complex JC(Bq) is bounded by k 2 ° { ' . 

Proof. The proposition is a consequence of Proposition [5] and the fact that the 
number of cells in the complex C(Ag) is bounded by k 2 ° {l . □ 

We now define, 

M Q = C'(K(B Q ), 

where C'{K.{Bq) is the cellular co-chain complex of the regular cell complex K,(Bq). 

Let T~1(Tq) (resp. 'H(Bq)) be a suitably fine cellular subdivision of Tq (resp. 
Bq) and let 

&.Q : C(H(B Q )) ^ C.(H(Tq)), 
be the homomorphism induced by a cellular map, which is a cellular approximation 

of (j>2, Q- 

Let 4>q : |/C(-Bq)| — > Bq denote the homotopy equivalence shown to exist by 
Proposition [6] above and let 

4>' Q :C.{K'{B Q ))^C.{U{B Q )), 

be the homomorphism induced by a cellular approximation to where KJ (Bq) 
is a cellular subdivision of the complex IC(Bq). 

Since, each cell of K,(Bq) is a union of cells of IC'(Bq), there is a natural homo- 
morphism 

Q :C.(K(B Q ))->C.{lC'{B a )) 
obtained by sending each p-dimcnsional cell of IC(Bq) to the sum of p-dimcnsional 
cells of JC'(Bq) contained in it, for every p > 0. It is a standard fact that 9q and 
its dual, 0q, are quasi-isomorphisms. 
Let 

4>q = 6 q o fi Q o 0' 2 Q : C'(H(T Q )) -» C(K(B Q )), 
where 4>'q (resp. <fi' 2 q) is the dual homomorphism of 4>'q (resp. <p' 2 q). 
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Proposition 8. For < i < k — 1, the induced homomorphisms, 
r Q : H\C'{H(T Q ))) ^ H\M' Q ) 

are isomorphisms. 

Proof. The proof is clear since ipQ is a composition of quasi-isomorphisms. □ 
Now let, B cAcV with #A = #B + K fc. 

The simplicial complex Ag is a subcomplex of A A and hence, /C(-Bg) is a sub- 
complex of fC{B A ) and thus there exists a natural homomorphism (induced by 
restriction), 

<t>A,B ■ M' A M* B . 
The complexes A4 A ,A4g, and the homomorphisms, 4>a,B,'4 ! A,4'B satisfy 



Proposition 9. The diagram 



... <PA,B ... 
M' A - 



■0^ 



IpB 



C'{H{T A )) * C'(H(T B )) 

is commutative, where r is the restriction homomorphism. 

Proof. Clear from the construction. □ 

It follows from Proposition [5] that the diagram (|6.ip is also commutative. 
We denote by 

fo,A ■ M' B -> M' A 

the homomorphism dual to <f> A g. We denote by Vfe* the double complex defined 
by: 

2>™ = M q 



The vertical differentials, 



Q- 

QCV,#Q=p+l 



are induced componentwise from the differentials of the individual complexes -Mg. 
The horizontal differentials, 

are defined as follows: for a <E = ©#g =p +i.Mg, and for each subset 

Q = {P io ,...,P ip+1 }cV 

with io < ■ • • < ip+i, the Q-th component of 6a £ Rp +1,9 is given by, 

0<J<P+1 
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1" 




i" 


— > 


I« 






— ► 


© #s =i.m! 




©#S=2.M 

1'. 


— > 


®#Q=iMq 


6 


©#S=2.M 


— > 


© #Q =iA ; 1q 



6 


©#s=2.M 




where Q i = Q\{P i .}. 



d 

e ©#q=3.Mq — ► ■•■ 

1" 

1' 



We have the following theorem. 

Theorem 6.2. For < i < k, 

IP(S) H\Tof(V^')). 

Proof. For each i, 1 < z < s, let C S k denote, the set defined on S k by P t < 0. 
Then, 5 = nf =1 Si. Choosing a suitably fine triangulation of U| =1 5i consider the 
Mayer- Vietoris double complex, A/"*'*, as described in Section I3~5l The homomor- 
phisms, 

-M 9 S ^ C q (H(T Q )) 

QCV,#Q=p+l Q C V,#Q=p+l Q C V,#Q=p+l 

give a homomorphism of the double complexes, 
(6.2) ^:V' V '^M'''. 

By Proposition [51 ip induces isomporphisms between the E\ terms of the two spec- 
tral sequences, obtained by taking homology with respect to the vertical differ- 
entials. Theorem 13.11 then implies that ip induces isomorphisms in the associated 
spectral sequences. But, the second spectral sequence converges to the homology 
of S by (|3.1| . The theorem is an immediate consequence. □ 

7. Algorithms for quadratic forms 

In this section we describe the algorithm for computing the top Bctti numbers of 
a basic semi-algebraic set defined by quadratic forms. We first describe an algorithm 
for computing index invariant triangulations (see Definition I6.1|) . 

Algorithm 1 (Index Invariant Triangulation). 

Input A set Q = {Pi, . . . , Pi } C R[A , . . . , Xk] where each Pi is a quadratic form. 
We denote by Pq = (Pi, . . . , P() : R fc+1 — > R £ the corresponding quadratic map. 
Output 

(A) A semi-algebraic triangulation, 

hn ■ An — > fto 
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of fig, which is compatible with the subsets f2g/ for every Q' C Q, and 
such that for any simplex a of Ag, index(a;P), as well as the multiplicities 
of the eigenvalues of wP, stay invariant as to varies over /ig(d). 
(B) For each simplex a of Ag, a set of parametrized univariate representa- 
tions, {uq(<j, lu), . . . ,Ufc((T, w)}, parametrized by u> G /ig(c), such that the 
associated points in R fc+1 , 

{e (<7,w), . . . ,e k (a,uj)}, 

form an orthornormal basis of R fc+1 , and 

{e 4 (cr,w), . . . ,e k (a,uj)} 

span the linear subspace of R fc+1 on which the quadratic form uiPq is pos- 
itive semi-definite (here i = index(wPg)). 
Procedure 



Step 1 Let e > be an infinitesimal and let Z = [Z\, . . . , Zi). Also, let Mq be the 
symmetric matrix corresponding to the quadratic form (in Xq, . . . , X k ) defined by 

P Q (Z) = (1 - e){ZiPx + ■■■ + Z e P e ) + sQ, 

where Q = Y)j—Q iXf. The entries of Mq depend linearly on Z\, . . . , Zi, and e. 
Compute the polynomial, 

F{Z, T) = det(M Q + T ■ h+i) = T k+1 + C k T k + ■ ■ ■ + C , 

where each d G R[e][Zi, . . . , Zg\ is a polynomial of degree at most k + 1. 

Step 2 Using Algorithm 11.1 in [9] (Elimination), compute a family of polynomials 
Aq C R[Zi, . . . , Ze] such that such that for each p G Sign(^lg), and z G H{p), 
the Thorn encodings of the roots of F{z 1 T) in R(e), as well as the the number of 
non-negative roots of F{z, T) stay constant. 

Step 3 Using the algorithm implicit in Theorem l3.2l fTriangulation). compute a semi- 
algebraic triangulation, 

: Ag — > Og, 
respecting the family Aq U L)j =1 {Zi}. 
Step 4 For each simplex a of Ag, let to(ct), . . . , Tfc(tr) be the Thorn encodings of the 
real roots, \o(cr,z) < ••• < Xk(cr,z) of F(z,T), for z G ^g(c). For < i < k, 
compute using linear algebra a parametrized univariate representation 

Ui(Z,T) = {F,g it0 , ■ ■ .,gi,k+i), 

such that for each z G ft-g(o") the point associated to Ui(z,T) is the eigenvector of 
Mq(z) corresponding to the eigenvalue Xi(z). Let u 2 ; = lim e Ui, and let Si(z) G R fe+1 
denote the corresponding unit vector. 

Complexity Analysis: The complexity of the algorithm is dominated by the 
complexity of Step 3, which is fc 2 1 '. □ 

PROOF of Correctness: The quadratic form Pq(z), and hence the matrix 
Mq, has k + 1 distinct eigcn values for each z G R f and has the same index as 
zPq. To prove the first part, replace e in the definition of Pq(Z) and observe that 
the statement is true t = 1, since the quadratic form Q has distinct eigenvalues. 
Thus, the set of t's for which Pq(z) has k + 1 distinct eigenvalues is non-empty, 
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constructible and contains a open subset, since the condition of having distinct 
eigenvalues is a stable condition. Thus, there exists £o > 0, such that for all 
t £ (0, £0)1 Pq( z ) nas k + 1 distinct eigen- values, and hence it is also the case for 
any infinitesimal t. The fact that Pq{z) has the same index as zPq follows from 
the fact that the quadratic form Q is positive semi-definite. Thus, the eigenspaces 
corresponding to the eigenvalues of Mq(z) are all one-dimensional and thus the 
vectors e^(er, u>) all well defined (upto multiplication by —1). Finally, note that 
since eg(er, lo), . . . ,e' k (a, to) are orthonormal, so are eo(cr, o>), . . . , efe(cr, u>) for every 
u) £ /ig(c), and since e[{a, ui), . . . , e' k (a, uj) span the non-negative eigenspace of 
luPu,, their images under the lim e map will span the non-negative eigenspace of 
loPq. □ 

We now describe an algorithm for computing the complexes A4q described in 
the previous section. 

Algorithm 2 (Build Complex for Unions) . 

Input (A) An integer £, < £ < k. 

(B) A quadratic map P = (Pi, . . . , P s ) : R fe+1 — > R s given by s homogeneous 
quadratic polynomials, Pi, . . . ,P S £ R[Xo, . . . , X^\. 
Output (A) For each subset Q C V = {Pi, . . . , P s }, #Q < £ + 2 a description 
of the complex M.q, consisting of a basis for each term of the complex and 
matrices (in this basis) for the differentials. 
(B) For each Q' C Q, with = j^Q! + 1, matrices for the homomorphisms, 

4>Q',Q ■ M' Q , -> M' Q . 

Procedure 

Step 1 For each subset Q = {P h P ie+2 } CP, with #Q = £ + 2, let P Q be the 
quadratic map corresponding to the subset Q. Call Algorithm Qjlndex Invariant 
Triangulation) with input Q. 

Step 2 Construct the cell complex C(Aq) (following its definition given in Section 

Step 3 For each cell C £ C(Ag), compute using the algorithm implicit in Theorem 
13.21 (Triangulation). the cell complex JC{C) and thus obtain a description of IC(Bq). 

Step 4 Compute the matrices corresponding to the differentials in the complex A4q = 
C'(IC(B Q )). 

Step 5 For Q' C Q C V with #Q = #Q' + 1<A:, compute the matrices for the 
homomorphisms of complexes, 

4>Q',Q ■ M' Q , -» M' Q . 

in the following way. 

The simplicial complex JC(Bqi) is a subcomplex of K,(Bq) by construction. Com- 
pute the matrix for the restriction homomorphism, 

e ,c '.C{K{B q ))^C{K{Bq,)). 

and output the matrix for the dual homomorphism. 

Complexity Analysis: The complexity of Step 1 YJiS> G)^ ^' usm g tnc com - 
plexity of Algorithm [TJ The complexity of Step 2 is Yli=o (i)^ 2 ° ( ' usm S the 

complexity of the algorithm for triangulating semi-algebraic sets. It follows from 
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Proposition [7] that the complexities of all the remaining steps are also bounded by 

Proof of Correctness: The correctness of the algorithm is a consequence of 
the correctness of Algorithm [1] and Propositions [6] and □ 

Let P±,...,P S G H[Xq, . . . , Xk] be homogeneous quadratic polynomials, and 
consider the set S C S k defined by S = {x € S k \ P 1 (x) < 0, . . . , P s < 0}. 

We will also denote for 1 < i < s, by Si the set defined by {x G S k \ P 2 (x) < 0.} 

Clearly S = f~\i<i< s Si. 

Algorithm 3 (Computing the highest I Betti Numbers: the homogeneous case). 

Input A quadratic map P = (Pi,...,P s ) : R k+1 — > R s given by a set, V = 

{Pi, ■ ■ ■ , P s } C R[Ao, . . . , Afc], of s homogeneous quadratic polynomials. 
Output b k (S), . . . , bk-e(S), where S is the set defined by 

S = p| {x g S fe | P(x) < 0}. 
Pep 

Procedure 

Step 1 Using Algorithm [2] compute the truncated complex £>*'*, i.e. 

2^,9 = 0<p<^+l, fc-i-l<g<fc, 

= 0, otherwise, 

Step 2 Compute using linear algebra, the ranks of 

JP(Tbt* (££'")), fc-*+l<i<&. 

Step 3 For each i, k-£<i<k, output, ^(S) = rank(.ff i (Tot*(X>]'"))). 

Complexity Analysis: The number of algebraic operations is clearly bounded 
by X)f=o (i)^ 2 ° ( * " usm S the complexity analysis of Algorithm [2l □ 

Proof of Correctness: The correctness of the algorithm is a consequence of 
the correctness of Algorithm [2] and Theorem 16.21 □ 

Remark 7.1. Suppose that (using Notation from Algorithm [3]) V' C V and 

S' = p| {x G S fe | P(x) < 0}, 
Pep' 

and letting T>''* denote the corresponding complex for S', it is clear from the 
definition that there is a homomorphism, Q-p^-pi : T>*'* — > ~D*'* defined as follows. 
For 

QCV,#Q=p+l Q C P,#Q=p+l 

§T>,-p-{(t>) = <fa- 

QCP',#S=P+1 

Recall from 16.21 that there exists, 

i> ■. v\' — ► a/-;** 

which induces an isomorphism, ^* : ff*(Tot*(2?*"*)) — ► H* (Tot' (J\f ''')). 
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Denoting by M'*'* the (truncated) Mayer- Vietoris complex for S' and by i-p,v '■ 
Afp'* — > Afp*'* the inclusion homomorphism, we have the following commutative 
diagram. 



H*(Tof(V'>')) 



V.V' 



H*(Tot'(V''')) 



H*(Tot'(Af;-'))) H*{Tot'(Af''''))) 

Note that H* (Tot' (Af*'')) H»(S) and F*(Tbt*(JV'*'*)) = H*(S'). 

It is clear that Algorithm [3] can be easily modified to output the complex £>*'*, 
by outputting the matrices corresponding to the vertical and horizontal homomor- 
phisms in the chosen bases. Furthermore, given a subset V' C V , Algorithm [3] 
can be made to output both the complexes and 2?/'* along with the matrices 
defining the homomorphism Q-p,v with the same complexity bounds. 

8. The General Case 

Let V = {Pi, ■ . ■ , PJ C R[Xi, X k ] with deg{P t ) < 2, f < i < s, and let 
S C R fe be the basic semi-algebraic set defined by Pi < 0, . . . , P s < 0. 

Let e > be an infinitesimal, and let P s+ i = £^j=i ~ 1- Let S C R(e) k be 
the basic semi-algebraic set defined by, Pi < 0, . . . , P a < 0, P s +i < 0. 

Proposition 10. The homology groups of S and S are isomorphic. 

Proof. This is a consequence of the conical structure at infinity of semi-algebraic 
sets (see for instance Proposition 5.50 in [9j)y. □ 

Moreover, denoting by P/ 1 the homogenization of Pi, and S h C S k the set defined 
by the system of quadratic inequalities, 

P[ l < 0,..., P s h <0,P s h +1 <0, 

on the unit sphere in R(e) k+1 we have, 

Proposition 11. For < i < k, b t (S) = \b l (S h ). 

Proof. First observe that S is bounded, and S h is the projection from the origin of 
the set S c {1} x R(e) k onto the unit sphere in R(e) k+1 . Since, S is bounded, the 
projection does not intersect the equator and consists of two disconnected copies 
in the upper and lower hemispheres, and each copy is homeomorphic to S. □ 

Algorithm 4 (Computing the top I Betti Numbers: the general case). 

Input A family of polynomials {Pi, . . . , P s } C R[Xi . . . , Xk], with deg(P^) < 2. 
Output b k _i(S), . . . , b k _i(S), where S is the set defined by 

S = f]{xe R k | P(x) < 0}. 
Per 

Procedure 

Step 1 Replace the family V by the family, V h = {Pf, . . . , P s '\ P s ' l + i}- 
Step 2 Using Algorithm [3] compute b k (S h ), . . . , b k -i(S h ). 
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Step 3 Output b k (S) = ±b k (S h ), b k ^ e (S) = \b k ^{S h ). 

Proof of Correctness: The correctness of Algorithm [4] is a consequence of 
Propositions [10] and [Tl] and the correctness of Algorithm [3] □ 

Complexity Analysis: The complexity of the algorithm is clearly 

i+2 

E 

from the complexity analysis of Algorithm [3] □ 

As in the case of Algorithm [3jsee Remark 17. 1[) Algorithm [4] can be modified 
to output a complex whose associated total complex is quasi-isomorphic to the 
truncated Mayer- Vietoris complex of S. Furthermore, given a subset V' C V 
defining the set S' D S, Algorithm [4] can be made to output the corresponding 
complexes of both sets, S and S' , as well as the homomorphism between them 
induced by inclusion. Since this fact is important in certain applications we record 
it as a separate theorem. 

Theorem 8.1. There exists an algorithm, which takes as input a family of poly- 
nomials {Pi, ■ ■ ■ , P s } C R[Ai . . . , Xk], with deg(Pi) < 2, and a number I < k, and 
outputs a complex T>*'* . The complex Tot*(2?* : *) is quasi-isomorphic to C,(S), the 

truncated singular chain complex of S, where 

S = p| {x g R k I P(x) < 0}. 
Per 

Moreover, given a subset T" C V , with 

S' = p| {x G R fe I P(x) < 0}. 
Pev 

the algorithm outputs both complexes and T>''' (corresponding to the sets S 
and S' respectively) along with the matrices defining a homomorphism Q-p.v, such 
that <5>* v v , : #,(Tbf (£>]'•)) = H*(S) H*(S') = i?*(Tot'(D'*'*)) is the homor- 
morphism induced by the inclusion i : S > S' . The complexity of the algorithm is 



9. Hardness 

In this section we show that the problem of computing the first few Bctti numbers 
of a semi-algebraic set defined by quadratic inequalities is #P-hard. Note that 
PSPACE-hardness of the problem of counting the number of connected components 
for general semi-algebraic sets were known before |T3l [26] and the proofs of these 
results extend easily to the quadratic case. 

Recall that in the classical Turing machine model, a function taking its values 
in N belongs to the class #P, if there exists a polynomial time non-deterministic 
Turing machine, whose number of accepting paths on any particular input is equal 
to the value of the function on that input [25] . A function / is #P-hard if every 
function in ffP can be reduced in polynomial time to the computation of /. The 
problem of counting the number of satisfying assignments to a Boolean formula is 
a #P-hard problem [25] . 
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The results of this section complements the main algorithmic result of this paper 
proved in the last section, that the problem of computing the top few Betti numbers 
of such a set is in P. 

Theorem 9.1. Given a family of polynomials V = {Pi,...,P s } C R fc , with 
deg(P) < 2,P G V, as input, the problem of computing be(S), where S C R fc 
is defined by P x > 0, . . . , P s > and I = 0(log£;) 7 is #P-hard. 

Proof. We first prove the case when I = by a straightforward reduction from 
the Boolean satisfiability problem. Clearly, if {C±, . . . ,C m } is an instance of the 
Boolean satisfiability problem with clauses C\, . . . , C m and variables Z±, . . . , Z n , 
then the number of satisfying assignments is equal to the number of connected 
components of the set defined by, 

Xi(Xi-l) > 0,Xi(l -XJ >0,...,X n (X n -l) >0,X n (l-X n ) >0 
C x - 1 > 0, . . . , C m - 1 > 0, 

where C; £ H[Xi, . . . , X n ] is the linear polynomial obtained from the clause Cj by 
syntactically substituting all the disjunctions by additions, and for each j, 1 < j < 
n, the literal Zj by the variable Xj, and the literal ~^Zj by the expression (1 — Xj). 

For I > 0, we reduce to the case £ = using the following observation. Given 
a basic semi-algebraic set S C R fe defined by, Pi > 0, . . . , P s > 0, deg(Pi) < 2, 1 < 
i < s, one can define another basic semi- algebraic set S' C R N defined by M 
polynomials inequalities with degrees < 2, with M = s+4m+l and N = k + 2m + l, 
where m is the total number of monomials in the polynomials Pi , . . . , P s . The 
number of monomials in the new system is bounded by 10to. The semi-algebraic 
set S' has the homotopy type of the suspension SS*. Moreover, the description of 
S' can be computed in polynomial time from the description of S. It follows from 
the basic properties of the suspensions (see [35]) that 6i(SS') = bo(S), which proves 
that computing bi(S) is also =#=P-hard. Iterating the construction, that is taking 
suspensions of suspensions £ times, and noting that be(Y< £ (S)) = bo(S), gives the 
result for £ = 0(log k) since the number of inequalities and variables increases only 
polynomially in s and k for £ = O(logfc). 

We now describe the construction of the set S". Introduce one new variable, Xq 
and consider the semi- algebraic set, 5"' C R k+1 defined by, 

(1 - Xl)Px > 0, . . . , (1 - X%)P S > 0, 
1 - xl > 0. 

The set S" = ([—1, 1] x S)L)HiUH2 where Pi and H2 are the hyperplancs defined by 
Xq = 1 and Xq = — 1 respectively It is easy to see that S" has the same homotopy 
type as the suspension of S. However, the polynomials used in the description of S" 
can have degrees as large as 4. We show below that by introducing new variables 
any quartic polynomial inequality can be reduced to a set of quadratic inequalities. 
The set S' is defined as follows. Write each monomial m of degree > 2 appearing in 
a polynomial used in the definition of 5", as a product of two monomials, mi, 1712, 
each of degree at most 2. Now replace each ocurrence of m in the inequalities used 
in the definition of S" , by the quadratic monomial Y m ,-\Y m ,i, where Y m ,\^ra,i are 
new variables. Also, add the inequalities 

Y m ,\ - m\ > 0, m x - Y m s > 0, 
Y m , 2 - m 2 > 0, m 2 - F m , 2 > 0, 
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to the set of inequalities. Clearly, the number of monomials, the number of new 
variables and the number of additional inequalities in the new system is bounded 
by 10m, 2m + 1 and 4m + 1 respectively, where m is the number of monomials in the 
original system. Finally, it is clear that the set S' defined by the above inequalities 
is a basic closed semi-algebraic set defined by polynomials inequalities of degree at 
most 2, which is homcomorphic to S", and hence has the same homotopy type as 
S5. This completes the proof. □ 
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